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ABSTRACT 


The parameter estimation of the Marshall- Olkin extended Burr III (MOEBIII) distribution, 
which is a generalization of the Burr III distribution is considered. The maximum likelihood (ML) estimation of 
the parameters of the MOEBIIIdistribution is introduced by Al-Saiarief al. (2016). However, ML is often used to 

estimate the parameters of the Burr III distribution; this method is very sensitive to the presence of outliers in the data. 

This paper presents M-estimation as a robust method, based on the quantile function to estimate the parameters of 
the MOEBIII distribution for complete data with outliers. A simulation study and areal data are used to illustratethe 
performance of M-estimator and ML estimator. The numerical results show that the M-estimator, generally is appropriate 
than the ML, in terms of the bias and mean square error when there are outliers in the data. 

KEYWORDS :Marshall- OlkinExtended Burr III, Maximum Likelihood, M-estimator, Robust estimator, Outliers 

1. INTRODUCTION 


In(1942) Burr introduced a family of distributions, in which, the Burr type III distribution is one of this family. 
The Burr type III distribution covers a wider region in the skewness and kurtosis plane.including several distributions such 
as the gamma family, Weibull family, lognormal family, and the bell-shaped beta distributions, see Tadikamalla (1980). 
It has been widely used as a model in many areas such as environmental, economics and reliability analysis. 

Marshall and Olkinin (1997) proposed a new family of survival functions based on adding a parameter to obtain a 
family of distributions, this distribution is called Marshall -Olkin extended distribution. 

Using the survival function of any distribution F(x), the survival function of the new distribution G(x; a) can be 
defined as follows: 


G(x] a ) 


a Fix ) 

1-a.Fix) 


— OO <X<CO,0<a<CO 


( 1 ) 


wherea = 1 — ais an additional positive parameter and they have called the tilt parameter 
[See Marshall and Olkin (2007)]. The cumulative distribution (cdf) and the probability density function (pdf) for the new 
distribution are given respectively by: 


G(x; a ) 


fQc) _ _ aFfr) 

1 -aF(x) l-aF(xy 


(2) 


g(x,a) 


a fix) 

[1 -aFix)] 2 


(3) 
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There have been various studies reported in the literature dealing with the parameter estimation methods for Burr 
III distribution. Lindsay et al. (1996) compared between the Burr type III distribution with four parameters and the Weibull 
distribution using the diameter data. Shao (2000) investigated the constrained ML method to estimate the parameters of 
the Burr type III distribution for toxicity data. Shao et al. (2008) used of the extended Burr type III distribution under 
the re-parameterization method in low-flow frequency analysis, where the lower tail of the distribution is of interest and 
used three methods of estimation, method of moments, probability-weighted moments, and maximum likelihood method to 
estimate the unknown parameters of the distribution. 

An alternative robust estimation method based on M-estimators and AM-estimators estimation method for 
the parameters of Burr III distribution have been proposed by Wang and Wen Lee (2010, 2011 and 2014) respectively as 
follows: 

In (2010) theyused the M- estimator to estimate the parameters of the extended three- parameter Burr type III 
distribution for complete data with outliers and compared it with ML and least squares methods. In (2011) they introduced 
M- estimator and ML and AM- estimator, when the data is asymmetric to estimate the parameters of the Burr type III 
distribution for complete data with outliers and compared the results. Also, in (2014) they discussed the M-estimator for 
estimating the Burr type III parameters with outliers and compared it with ML and least squares methods. 

However, there is not much work concerning MOEBIII distribution. The parameters of MOEBIII distribution was 
estimatedby using ML estimation method by Al-Saiariet al. (2016). However, it is well established that, in the presence of 
outliers in the data, the traditional methods do not provide reliableestimators. Therefore, robust estimation methods can be 
used for estimating the parameters of the MOEBIII distribution, if the data contains outliers. 

This paperis organized as follows: In Section 2, the Marshall- Olkin extended Burr III (MOEBIII) distribution is 
described.In Section 3, the maximum likelihood and therobust M-estimator are introduced. A simulation study and a real 
data to illustrate the performance ofML estimator and M-estimator are illustrated in Sections 4 and5, respectively. Finally, 
concluding remarks are given in Section 6. 

2. MARSHALL-OLKIN EXTENDED BURR III DISTRIBUTION 


The cumulative distribution (cdf) and the probability density function (pdf) for a two shape parameters c > 0, 
k > 0 Burr type III distribution are given by: 

Fix; c, k) = (1 + x~ c T k ; x > 0 (4) 

and 


fix; c, k) = kcx< c+1 \ 1 + x- c )- (k+1) ; *>0 
where c, k are shape parameters. 

Substituting (4) and (5) in (2) and (3) the cdf and the pdf of MOEBIIIare given respectivelyas: 

(: i+x~ c )- k 


Gix; a, c, k ) 


[l-(l- a )[l-(l+*- c r fc ]] 


; x > 0 


(5) 


(6) 
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and 


g{x ; a, c, k ) 


akcx~( c+ 1 \l+x~ c )~( k+1 ) 

2 ; X > 0 

[l-(l- a )[l-(l+x- c )- k ]] 


(7) 


where, a is called tilt parameter, whena = 1, the original distribution i.e., Burr type III distribution with two 
parameters is obtained. 


MOEBIII distribution is more flexible than the BurrIII distribution, because of the presence of the additional 
shape parameter. Figure 1 shows the plots of pdf for MOEBIII distribution for some values of the parameters 
[See Al-Saiari .et.al. (2016)] 



Figure 1: Plot of the PDF of the MOEBIII distribution. (I) a = 0.8,k=0.1,c=2, 
(II) a = 3, k=0.1,c=2, (III) a = 0.8, k = 5, c = 2, (IV) a = 3, k = 5, c= 2 


The qthquantile of the MOEBIII is given by 


x q = G 1 (q) 


[fir) + 1 


v k 


; 0 < q < 1 


3. PARAMETERS ESTIMATION 


(8) 


In this section, MLestimation method to estimate the parameters of theMOEBIII distributionis first 
reviewed. Therobust estimators based on the M-estimation method given by Huber (1964) are introduced. 

3.1 The Maximum Likelihood Method 

Suppose thatA = (x 1 ,x 2 , ... ,x n ) is a random sample of size n from MOEBIII(a, c, k), then the log-likelihood 
function I (a, c, k ) is given by: 

l(a, c, k ) = n(loga + logc + logk) — (c + 1) £"= t log(Xi) — (k + 1) J]f =1 log( 1 + x 1 rC ) — 2 £f=i 1 — 

(1 - a) [1 - (1 + xr c r k ] = 0 (9) 
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To obtain ML estimates of the parameters, equation (9) is maximized by taking the partial derivatives with respect 
toa, candk setting it to 0 as follows: 


3l{a,c,k ) _ n n v , n 8 ) _ n 

■ - » z ~~ u 


da cL 


i-(i-a)[i-(i+*r £ ) fc ] 

dl(a,c,k ) n v -1 " 


dc 


Tl \~~* n y ^ 71 

- - / logixi ) + (fc + l) > 

c £ — H =1 


n X t e log(x l ) 
=i (1 + x ; - e ) 


2(1 - a)k’Zf =1 Xi ~ Cl09iXi ^ 1+Xi ~ C) [ J = 0 


i-(i-2)[i-(i+^r 8 ) k ] 


( 10 ) 


(ii) 


d(a,c,fc) n 


dk 


^-Sf =1 /o 5 (l+xr e ) + 2(l-a)Ef =1 


n ( 1+ *i c ) l °a( 1+Xj 8 ) _ 


i-Ci-a)|i-(i+Ar i £ ) k j 


( 12 ) 


Since equations (10) - (12) are non-linear equations, it cannot be solved analytically;numerical methods 
shouldbeused to obtain the ML estimates. 


3.2 M- Estimation Method 

The robust M-estimator methodis used to estimate the parameters of the MOEBIII distribution from a sample of 
nobservations(x 1 ,x 2i ... , x n ) = x T . As follows: 

x i = /;(#) + u t ;i = 1, ... ,n (13) 

where# 7 = ( a , c, k ), /;(#) is the quantile function of the MOEBIII distribution andtt; is the error term with mean 
equals zero and varianceequalser 2 . Gilchrist (2000) introduced the robust property based on the quantilefunction. 
For the robust regression method, it is necessary to scale the invariant error, as follows: 



wheres = and MAD (it) = Mi4D(u 1 ,it 2 , ...,ti n ) = Median{\u — Medianiu) |) 

[See Hampelet.a/. (1986)] 

The M-estimator method for estimating the MOEBIII parameters is defined by minimizing the objective function 
of all invariant errors, p(e ; )as follows: 

Minimize £f=iP( e i)(15) 

To estimate the three unknown parameters of the MOEBIII distribution, a simple comparison among two different 
p objective functions is used. The selected objective functions areTukey’sBisquare and Huber’s weight 
[See Huber (1981)]. 

Tukey’sBisquare objective function is derived as: 
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Pied = 


, JtM-(3 2 )Vw 


< a 


-a 2 if\e i \ > a 


with the derivative 


’/'©;) = { e i 


(l-© 2 ) lf\e,\ 


0 


< a 
if |e ; | > a 


where a= 4.685 ande;is the scale invariant error in (14). 
Huber’s weight objective function is derived as: 


P( e i) 


-(e;) z t/k| < a 
a |e £ | — -a 2 if\e i \ > a 


with the derivative 

Me ) = tf\ e i\ < a 

^ l} 1 asigne t if \e t \ > a 

where a= 1.345 ande ; is the scale invariant error in (14). 


( 16 ) 


(17) 


The estimators of the parameters can be obtainedfor the two selected objective functions are derived by 
differentiating (15) with respect toa, cand/crespectively, and equating to zero. Then, the simultaneous equations can be 
obtained, as follows: 


M 


-1 


r/1 - q\ 

i 

k 

+ 1 

- 1 

L\ aq ) 



-(H 


[© + f '] = 




+ 1 


- 1 


x log 


[©) 


+ 1 - 1=0 




(— ) 

V aq ) 


+ 1 


- 1 


- 6*0 


((©) +1 ]'M(©) +1 ] = 0 


(19) 


( 20 ) 
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Since equations (18)-(20) are non-linear equations, numerical iterative methods can be used to obtain the M- 
estimators for the parameters a, candk. 

4. SIMULATION STUDY 

In this section, a simulation study to compare the performance of the MLand the robust M-estimatorsis carried out 
in the presence of outliers. The data were generated from the MOEBIII distribution for several different values of 

9 T =a, c and k. 

The steps of simulation study are given as: 

• Step 1. Generate n random samples w 1 ,w 2 , ... , w n from uniform distribution ~U (0, 1). 

• Step 2 Generate a random sample*; from the MOEBIII distribution for specified values of a, candk distribution 
using equation (8) for different valuesoft, i = 1,2, ... ,n. 

• Step 3. Select different sample sizes, n= 20, 40,100 to represent the small, moderate and the large sample sizes, 
respectively. 

• Step 4. The outliers are generated for each random sample from the uniform distribution as U (X+ 4S, X + IS), 
whereXthesamples mean ofisX = (x 1 ,x 2 ...,* n )andS is the sample standard deviation of X. The outliers are 
generated by shifting the largest observations to the right in the X direction. For the small sample size (n = 20) 
one outlier is taken, for the moderate sample size (n = 40) two outliers are taken and for the large sample size (n = 
100) five outliers are taken. [See Wei and Fung (1999)]. 

• Step 5. Different initial parameter values for a, c and k are taken. 

• Step 6.0btain the ML estimates using equations (10) - (12). 

• Step 7,Obtain the M- estimates using equations ( 1 8)-(20). 

• Step 8. calculate the bias and mean square errors (MSE) of the estimates derived in the steps 6 and 7 using the 
following formulae: 

Bias{6 ) = E{G ) - 9 

N 

i= 1 

where, E(0) = 

• Step 9. Repeat steps (1-7), N=1000 time. 

The results in Tables 1-3 clearly indicates that, the robust estimator based on Tukey'sBisquarefunction has the 
smallest bias and the smallest MSE in most of the cases are with outliers. In addition, the robust estimator based on 
Huber’s function outperforms ML estimator in terms of bias and MSEwhen the data set contains outliers. For example, it 
can be seen from Tables 1-3 that, the largest difference of the bias for parameter c arises in Table 1 for the cases (0.2, 2.0, 
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0.2) and (0.3, 2.0, 0.1). Similarly, one can observe superiority of the robust estimators in terms of bias and MSE in Tables 
1-3. Moreover, in Table 3 it can be noted that when the number of outliers increase, the ML estimators increase when 
compared to the robust estimators. It can be seen that, bias and the MSE values are very large for the ML estimators. In 
general, it is advised to use the MOEBIII distribution instead of using classical methodswhen there are potential outliers in 
the data. 


Table 1: The Bias and MSE (in Parenthesis) for n=20 with One Outlier 


Parameter a 

Parameters 

ML 

Tukey’sBisquare 

Huber 

(0.1, 2.0, 0.1) 

0.2695 (0.1485) 

0.2735 (0.1057) 

0.2653 (0.1086) 

(0.2, 2.0, 0.2) 

0.2207 (0.1146) 

-0.0855 (0.0267) 

0.1357 (0.0732) 

(0.3, 2.0, 0.3) 

0.1097 (0.1038) 

-0.1758 (0.0419) 

-0.0239(0.0528) 

(0.3, 2.0, 0.1) 

0.2470 (0.1225) 

-0.0606 (0.0502) 

0.1884(0.0926) 

(0.8, 2.0, 0.1) 

-0.1333 (0.0742) 

-0.0279 (0.0043) 

-0.1647 (0.0682) 

Parameter c 


ML 

Tukey’sBisquare 

Huber 

(0.1, 2.0, 0.1) 

-0.3315 (6.2178) 

-0.47301 (0.7879) 

-0.0661 (0.2523) 

(0.2, 2.0, 0.2) 

1.3283 (41.1222) 

0.2037 (9.3444) 

-0.1857 (0.6362) 

(0.3, 2.0, 0.3) 

-0.2962 (5.4922) 

-0.3118 (1.5008) 

-0.62901 (1.0165) 

(0.3, 2.0, 0.1) 

1.1688 (32.7146) 

0.9089(11.9765) 

0.3605 (5.2050) 

(0.8, 2.0, 0.1) 

0.0501 (11.1791) 

-0.3109 (3.2648 ) 

-0.7966(0.8742) 

Parameter k 


ML 

Tukey’sBisquare 

Huber 

(0.1, 2.0, 0.1) 

0.6996(2.4519) 

-0.0316(0.0019) 

0.1019(0.2715) 

(0.2, 2.0, 0.2) 

0.0545(0.0638) 

-0.0772 (0.0122) 

0.0601(0.0151) 

(0.3, 2.0, 0.3) 

0.4316 (1.0632) 

0.0856 (0.2068) 

0.3009 (0.6877) 

(0.3, 2.0, 0.1) 

0.0418(0.0296) 

-0.0083(0.0086) 

0.0167 (0.0078) 

(0.8, 2.0, 0.1) 

0.0427 (0.0159) 

-0.0279 (0.0043) 

0.03561(0.0055) 


Table 2: The Bias and MSE (in Parenthesis) for n= 40 with 2 Outliers 


Parameter a 

parameters 

ML 

Tukey’sBisquare 

Huber 

(0.1, 2.0, 0.1) 

0.1162 (0.0319) 

0.0876 (0.0293) 

0.1421 (0.0218) 

(0.2,2, 0.2) 

0.3606(0.1803) 

0.0022 (0.0556) 

0.1483 (0.0615) 

(0.3,2,03) 

0.1445(0.0904) 

-0.0557 (0.0474) 

0.0327 (0.0527) 

(03,2.0,0.1) 

0.2459 (0.1024) 

-0.0031 (0.0651) 

0.1805(0.0729) 

(0.8, 2.0, 0.1) 

-0.4973 (0.3170) 

-0.1038 (0.0552) 

-0.3035(0.1660) 

Parameter c 
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ML 

Tukey’sBisquare 

Huber 

(0.1,2. 0,0.1) 

0.2669 (11.7336) 

-1.2805 (1.7296) 

-0.9759 (1.1687) 

(0.2,2, 0.2) 

-1.2267 (1.7687) 

-0.6979 (0.7243) 

-0.5429 (0.4339) 

(0.3,2,0.3) 

-0.3030 (1.2793) 

-0.9431 (1.1152) 

-0.8058 (1.2048) 

(03,2.0,0.1) 

0.3482(11.8764) 

-0.4349 (3.1586) 

-0.3435 (1.6942) 

Table 2: Contd., 

(0.8,2. 0,0.1) 

-0.4489(1.9466) 

-0.8765 (1.6520) 

-0.6734(1.1204) 

Parameter k 


ML 

Tukey’sBisquare 

Huber 

(0.1,2. 0,0.1) 

0.1368(0.0356) 

0.0386 (0.0083) 

0.1030(0.0202) 

(0.2,2, 0.2) 

0.10367 (0.0286) 

-0.0399 (0.0089) 

0.0076 (0.0027) 

(03,2,0.3) 

0.1186 (0.1028) 

-0.1806(0.0472) 

0.0555(0.0554) 

(03,2.0,0.1) 

0.0768 (0.0269) 

0.0162 (0.0208) 

0.0555 (0.0095) 

(0.8,2. 0,0.1) 

0.1057 (0.0435) 

-0.0074 (0.0054) 

0.0998 (0.0217) 


Table 3: The Bias and MSE (in Parenthesis) for n=100 with 5 Outliers 


Parameter a 

Parameters 

ML 

Tukey’sBisquare 

Huber 

(0.1, 2.0, 0.1) 

0.2608 (0.1478) 

-0.0199 (0.0004) 

0.1808(0.0696) 

(0.2, 2.0, 0.2) 

0.1126 (0.0530) 

-0.0453 (0.0375) 

0.0683(0.0239) 

(03,2.0,03) 

0.1207 (0.0581) 

-0.0991 (0.0319) 

0.04755 (0.0410) 

(03,2.0,0.1) 

0.2354 (0.0928) 

0.0067 (0.0442) 

0.1290 (0.0397) 

(0.8, 2.0, 0.1) 

-0.4048 (0.2387) 

-0.0191 (0.0205) 

-0.0249 (0.0213) 

Parameter c 


ML 

Tukey’sBisquare 

Huber 

(0.1, 2.0, 0.1) 

0.4732(4.2175) 

-0.0575 (1.9262) 

0.2640(1.50503) 

(0.2, 2.0, 0.2) 

-1.2290(1.8074) 

-0.7163 (0.5969) 

-0.7037 (0.8123) 

(03,2.0,03) 

-0.7874( 1.0200) 

-0.6787 (0.7803) 

-0.6569 (0.7127) 

(03,2.0,0.1) 

0.0947 (5.3904) 

-0.6648(1.1848) 

-0.2354 (2.1846) 

(0.8, 2.0, 0.1) 

-1.1227 (1.5925) 

-0.7149(0.6458) 

-0.7023 (0.8461) 

Parameter k 


ML 

Tukey’sBisquare 

Huber 

(0.1, 2.0, 0.1) 

0.7522 (3.1906) 

0.0112(0.0072) 

0.0776 (0.0258) 

(0.2, 2.0, 0.2) 

0.1934 (0.0923) 

0.0844 (0.0332) 

0.2331 (0.1117) 

(03,2.0,03) 

0.2987 (0.7993) 

0.0035 (0.2025) 

0.0965 (0.0278) 

(03,2.0,0.1) 

0.0621 (0.0219) 

-0.0179 (0.0033) 

0.0513 (0.0109) 

(0.8, 2.0, 0.1) 

0.0833 (0.0102) 

-0.0068 (0.0046) 

0.0606 (0.0089) 
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5. APPLICATION TO REAL DATA 


In this section, the data set given by Shao (2000) is used. This data set is used by Wang and Wen Lee (201 1 and 
2014), to illustrate the Burr type III distribution and the ability of the M- estimation method and AM- estimation method. 
This data are the result of the study of influence of the proportion of toxicityof chromium in marine waters. 


The data set, which is given in Table 4, contains 36 for chromium in marine waters. The box plot of the dataset 
given in Figure 2 displays a potential outlier in the data set. 


Table 4: Number Observed Effect Concentration (NOEC) Values for Chromium in Marine Water 


Marine Water 

2000, 12.59, 199.53, 4, 56, 263.03, 776.25, 1258.93, 3311.31, 39.81, 2630.27, 187.2, 89.13, 4.79, 8800, 10000, 
1600, 602.56, 2000, 199.5, 2.4, 728, 1122.02, 140, 177.83, 540, 1456, 2511.89, 1200, 478.63, 1174.89, 9.55, 3090.3, 
602.56, 1995.26,210. 



O 2000 4000 eooo SOOO lOOOO 


Figure2: Box plotof the Chromium in Marine Water Data 

The data set is used to fit a MOEBIII distribution and to estimate the unknown parameters a, c and k, using the 
ML and the M-estimator based on two different objective functions:Tukey’sBisquare and Huber’s function. Table 5 gives 
summary of fitting the MOEBIII distribution obtained from theML, and the robust estimation methods for this dataset. The 
95% confidence intervals of a, c andkfor ML and M-estimator with the two objective functions, Tukey’sBisquare and 
Huber’s are shown in Table 5. In Table 5, one can see that, the confidence intervals of Tukey’sBisquareare estimates are 
better than the other fitted obtained from the Huber and ML estimates.Also, it is found that Huber estimates is better 
thanthe ML estimates. 


Table 5: Parameter estimates, Standard Deviation (in Parentheses) and 
Confidence Intervals [in Parentheses] forthe data Set 

Parameters Method 
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ML 

Tukey’sBisquare 

Huber 

a 

5.1(0.852) 

[3.43,6.77] 

5.9(0.110) 

[5.784,6.215] 

5.6(0.273) 

[5.065,6.135] 

c 

0.65(0.0739) 

[0.40,0.795] 

0.71(0.035) 

[0.641,0.779] 

0.59( 0.0583) 
[0.476, 0.704] 

k 

7.8(1.705) 

[5.66,11.14] 

8.9(0.098) 

[8.708,9.092] 

8.3(0.197) 

[7.914,8.686] 


6. CONCLUSIONS 

The M-estimator method based on the quantile function by the Tukey’sBisquare and Huber objective function to 
estimate the MOEBIII parameters for complete data with outliers are presented. The simulation results showed that the M- 
estimator method outperforms the ML method regarding bias and mean square error. The real data showed that the M- 
estimator method is suitable for estimating the MOEBIII parameters for complete data with outliers. Also, the confidence 
intervals of the MOEBIII parameters for M-estimator and ML are obtained using this data set. The confidence intervals of 
Tukey'sBisquare estimates are comparably better than the other fitted, obtained from the Huber and ML estimates. The 
Huber estimates is also better than the ML estimates. In general, it is advised to use the MOEBIII distribution instead of 
using classical methods, when there are potential outliers in the data. 
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